Random Walks, Reaction-Diffusion, and Nonequilibrium Dynamics of Spin Chains in 

One-dimensional Random Environments 
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Sinai's model of diffusion in one-dimension with random local bias is studied by a real space 
renormalization group which yields asymptotically exact long time results. The distribution of the 
position of a particle and the probability of it not returning to the origin are obtained, as well as the 
two-time distribution which exhibits "aging" with scaling and a singularity at lnt = \nt'. The 
effects of a small uniform force are also studied. Extension to motion of many domain walls yields 
non-equilibrium time dependent correlations for the ID random field Ising model with Glauber 
dynamics and "persistence" exponents of fD reaction-diffusion models with random forces. 



The development of order in systems with a broken 
symmetry is of interest in many contexts. "Coarsening" 
of domain structures evolving towards equilibrium has 
been studied extensively But little is known analyt- 
ically about domain growth in the presence of quenched 
disorder Nevertheless, phenomenological descrip- 

tions of the non-equilibrium dynamics of various ran- 
dom magnetic systems have been developed in terms of 
" droplets" separated by domain walls || . Due to the very 
slow dynamics associated with activation over large free 
energy barriers, even the apparent equilibrium properties 
of these systems are dominated by the non-equilibrium 
dynamics, as also occurs in infinite-range models [pH. 
Even in one-dimension some random systems exhibit ul- 
tra slow growth and aging phenomena. Exact results in 
ID could thus be used as testing grounds for more com- 
plex D > 1 cases which have resisted analytic attack. 

In this Letter we study the diffusion of a particle in a 
ID random potential which itself has the statistics of a 
ID random walk Q. Extensions to many interacting par- 
ticles allows us to study, via domain walls, the Glauber 
dynamics of ID Ising models, in particular random field 
ferromagnets and spin glasses in a magnetic field. This 
leads also to the consideration of more general diffusion- 
reaction processes in such energy landscapes. Various 
analytic results are known for the single particle model 
(Sinai model) [f^-p"^ but these primarily concern single 
time quantities. Here we use a real space renormalization 
group (RSRG) method related to that used to study dis- 
ordered quantum spin chains H-M. This allows us to 
compute a host of quantities such as first passage (per- 
sistence) exponents, single time correlations and even 
two time correlations that are probed in aging exper- 
iments. Despite its approximate character, the RSRG 
yields many results that should be exact. 

The model is defined as follows: Particles diffuse on 
a ID lattice in a potential Ui, with i the site index. A 
"force" variable fi = Ui — L7 i+1 is defined on each bond 
+ 1) with these fi independent random variables. 
Since one can group together neighboring bonds with the 
same sign of the force, we study with no loss of generality, 



a " zigzag" potential (see Fig. [|) with the fi alternatively 
positive and negative but with a distribution of "bond" 
lengths k. Our model is thus defined by /, = (— l) l+1 Fi 
where the positive Fi = \Ui — Ui+i\, which are effectively 
energy barriers, are the natural variables. The pairs of 
bond variables F, I are chosen independently from bond 
to bond from a distribution P(F,l). In the presence of 
a directionality bias one needs two distinct distributions 
P(F, I) for "descending bonds" and R(F, I) for "ascend- 
ing bonds" , both normalized to unity. 
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FIG. 1. (a) Energy landscape in Sinai model (b) decima- 
tion method: the bond with the smallest barrier F m i n = F2 
is eliminated resulting in three bonds being grouped into one. 



We are interested in long times when the behavior will 
be dominated by large barriers and it is on these that 
we must focus. Our RG procedure is conceptually sim- 
ple: in a given energy landscape it consists of iterative 
decimation of the bond with the smallest barrier, say 
F-2 = U3 — U2 as illustrated in Fig. [l]. At time scales much 
longer than enp(F2/T), local equilibrium will be estab- 
lished between sites 2 and 3 and the rate for the walker 
to get from 4 to 1 will be essentially the same as it would 
be if sites 2 and 3 did not exist but 1 and 4 were instead 
connected by a bond with barrier F' = F\ — F 2 + F3 
and length I' = h + h + h- We thus carry out ex- 
actly this replacement which preserves the zigzag struc- 
ture and the larger scale extrema of the potential. With 



1 



r defined to be the smallest remaining barrier at a given 
stage of the RG, we eliminate the barriers in the range 
r < F < T + dT. The new variables are independent from 
bond to bond. Introducing the variable C = F — T one 
finds the following RG flow equations for the probabilities 
P(C, I, T) and R(C, I, T) (with < C < oo): 

(d r - d ( )P = R(0, .) *, P* C iP+ (P r - R%)P 

(dr - d c )R = P(0, .) *, R* cl R+ (iZj - P r )i? (I) 

Here *^ denotes a convolution with respect to £ only 
and with respect to both £ and I and we define 
P r = dlP(( = 0, 1, T) and similarly for R$. The 
dynamics implied by this RG is rather simple. Making 
the obvious identification of T = Thi(t/to) from Arrhe- 
nius dynamics, we see that at very long time scales the 
renormalized landscape consists entirely of deep valleys 
separated by high barriers. A good approximation to the 
long time dynamics is thus to place the walker at the bot- 
tom of the renormalized valley at scale T In t in which it 
was initially, since, with high probability, it will be near 
to that point |Q . Upon proper rescaling of space and time 
this approximation becomes exact as T tends to oo as was 
proven in ref. [Q for the unbiased case. It remains valid 
in the weakly biased case in the limit that the bias pa- 
rameter that controls the long time properties, /i, defined 
implicitly for the original model with unit length bonds 
by < exp(-fifi/T) >= 1, is very small (see (1010). 

The RG equations (|]) are identical to those derived for 
the random transverse field Ising chain (RTFIC) in |H| 
with the identification of In hk = p2fe and In Jk = P2fc+i 
in the RTFIC with the ascending and descending barriers 
respectively Jl^] . Thus duality in the RTFIC corresponds 
to reversing the average force. Criticality then corre- 
sponds to the zero drift case, while the Griffiths phase 
in the RTFIC |I5f| corresponds to the biased phase with 



zero velocity |ld| , |l if . The deviation from criticality pa- 
rameter S = ((m/i) — (In J)) /[var (In h) + var(lnJ)] 
is analogous at small 6 to fi/2. 

We consider first the long time dynamics of a single 
particle, starting with the symmetric, zero bias, case 
that has the same distribution for all the bonds; i.e., 
R = P. For large T , the distribution flows to one of a 
family of scaling solutions of the RG equations ([!]). The 
rescaled probability P(r),\) = {T 3 /a)P(r]T,XT 2 /a,T) in 
terms of the rescaled variables r\ = £/r, A = la/T 2 , 
when Laplace transformed in A to s, is found to be jl5[ ] 
P(rj,s) = (y/s/ sinh( v / s)) exp(— r]y/s coth( v / s)). The av- 
erage bond length is I — -^T 2 and we recover the scaling 
x ~ In 2 t M . The large scale variance of the potential 
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with L_a the distance from i 



to j is conserved by the RG, fixing a. 

The fact that the renormalized barrier distribution be- 
comes infinitely broad is the source of the exactness of 
our long time results. At late stages of the RG, the 
chances that two neighboring bonds have F's that are 



within order T of each other tends to zero for large T. 
Thus substantial errors that are introduced by assigning 
a particle to one of two almost-equal-depth neighboring 
valleys rather than splitting its distribution between the 
two valleys will occur rarely at long scales. Furthermore, 
any such error is wiped out by a later decimation which 
eliminates the two valleys in favor of a deeper valley. 
Since in a deep renormalized valley, the particle tends 
to be very close to the bottom on the scale of 7(r) |?],[|, 
we can obtain the scaled distribution of the position of a 
particle at time t that started at the origin at time zero, 
directly from P(l,T = Tint). We henceforth set T = 1 
and measure distances in units such that a = 1. The 
dynamics "within" one bond and errors made in early 
stages of the RG will generally only change the micro- 
scopic cutoff time in lnt's. 

The renormalized dynamics corresponds to moving the 
particle from its starting point (distributed uniformly on 
a bond) to the lower-potential end of the bond. The 
distribution of its position at time t averaged over the 
ensemble of random potentials is thus Prob(a:, i|0, 0) = 
]™,dlP(l,T)/l(T). With T = lnt, it takes the scaling 

form Prob(a;, i|0, 0) = i^jlii^i) where 
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With a reinserted, this coincides with Kesten's rigorous 
result H for a Brownian potential, as it should . 

But we can now generalize to the biased case with a 
small average potential drop per unit length 28 > 0. The 
RG flows Eq(|) now involve the two distributions R and 
P. The asymptotic behavior of these flows was found in 
IP^I ; it obtains in the scaling limit that F is large, while 
A = l/T 2 and 7 = TS are both fixed but arbitrary. In 
terms of the Laplace transform from \X\ = |^|/r 2 to s, 
we obtain for the generalization of Eq(||), 



9(X, 7 ) = (-^r-f[9(X)LT-'-(l 
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s k cosh k + 7 smn k 



k cosh k — 7 sinh k 
)}, (3) 



with k = \/s + 7 2 and the two terms arising from 
descending and ascending renormalized bonds, respec- 
tively. In the limit of small 7, the behavior reduces 
to the symmetric case Eq(||). But for large 7, i.e 
hit >> 1/5, the distribution is heavily concentrated to 
the right of the origin being simply Prob(x, i|0, 0) » 
6(x) exp(— x/x(t))/x(t) with the mean displacement 
x(t) fts t 2S /(4S 2 ), consistent with the small /i limit of the 
known |80 "Levy front" L^tj '(fi 2 x) 1/>Jj ), although the 
exponent \i of the anomalous diffusion, x ~ t^ is correct 
only to leading order in 6, due to corrections to scaling 
neglected in our RG. We find that the model renormal- 
izes onto a directed model with traps (ascending bonds) 
of "release time" distribution p(r) ~ t"( 1+ ^ O. 
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Our method also enables us to compute two time quan- 
tities, e.g. B(x,t,x' ,t') = Prob(xt|x't'|00) which con- 
tains information about the dynamics after the system 
has "aged" from t = to t', but the full calculation and 
result are too complicated to reproduce here. In the 
regime r = lnt and T' = hit' large with a = lnt/ lnt' a 
fixed number, we obtain a scaling form B ps B a (x, x') in 
the rescaled variables x = x/T 2 and x' = x' /T 2 . Our two 
time correlations thus exhibit a lnt'/ lnt aging regime, 
as found numerically in |l3| ] and argued in [|| for spin 
glasses in higher dimensions. Interestingly, the rescaled 
distribution has a delta function component at the origin, 
suggested in [|l3| but obtained here analytically. It arises 
from bonds not decimated between t' and t, i.e., from 
particles staying within the same valley. The probability 
Dr,r'(C) that a particle on a bond with F at V has not 
moved from V to T satisfies a linear RG equation which, 
when integrated, yields the probability that a walker has 
not moved substantially between t' and i 



(4) 



The mean square additional displacement (\x(t) — x(t')\ 2 ) 
at large a is ~ y^j(lnt) 4 . But when lnt and lnt' are not 
too separated, i.e a « 1 it is only (lnt') 4 |y|(lnt/ lnt' — 
1). Typically in this regime the particle is trapped in a 
deep well, but there is a probability of order (T — r')/r' 
that one of the barriers of the well at time t' is less than 
r. If this does occur, then the walker will jump to the 
bottom of a deeper valley a distance of order l(T') ~ T' 2 
away. Note, however, that for (t — t')/t' = 0(1) or less, 
this will not occur as t' — > oo and the behavior will in- 
stead be dominated by rare configurations -absent in the 
scaling limit- in which the valley at time t has two al- 
most degenerate minima. Jumping between such min- 
ima persists even for t — > oo with t — t' fixed and in 
this limit the statistics of the infinitely deep valley po- 
tential becomes that of a random walk restricted to have 
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We now turn to problems involving many walkers. The 
Glauber dynamics of the ID (classical) random field Ising 
model corresponds to two types of domain walls A and 
B which see opposite random potentials with the forces 
fi being simply twice the corresponding random fields 
on the dual lattice sites. When the random fields are 
much smaller than the exchange </, the long time behav- 
ior for T « J will be universal. We focus on the evolu- 
tion starting from random initial conditions - e.g., after 
a quench from a high temperature. At low temperatures, 
an A wall quickly falls to the bottom of a valley only to 
move to the bottom of a neighboring lower valley when 
lnt reaches the barrier height of the intervening bond. 
Likewise, B walls move from top to top of "mountains". 
When an A and a B meet, they annihilate preserving the 
alternating ABAB sequence. Analytic treatment is diffi- 
cult because decimation generates correlations among the 



A and B positions. Performing the RSRG numerically on 
a large sample [0 we find that the system evolves to a 
state with one A at each minimum and one B at each 
maximum of the renormalized landscape. We thus make 
the Ansatz that this is the correct form of the asymp- 
totic states. The RG analysis is then again simple Jl6| . 
The equal time spin correlations can be obtained from 
the difference between the probabilities that an even or 
an odd number of extrema of the renormalized potential 
-i.e., domain walls - exist between a given pair of points. 
For a symmetric distribution of random fields, we obtain: 
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with X = distances normalized as earlier and T = 
Tint. At sufficiently long times r > Tj — 2 J, we can 
no longer ignore creation of pairs of walls. But, at this 
point, the energy cannot be lowered further by any pro- 
cess. Thus if the renormalization is stopped at Tj, in the 
small field, low T scaling limit the configuration of the 
walls corresponds, up to negligible thermal fluctuations, 
to the equilibrium state. The above equation should then 
give the mean equilibrium spin correlation function with 
lengths measured in units of the Imry-Ma length above 
which the random fields dominate the exchange. 

Since ID Ising spin glasses in a field are equivalent via 
a gauge transformation to random field ferromagnets, we 
can also obtain results for such a system. If a large mag- 
netic field is quickly reduced to be << J but non-zero, 
the domain wall dynamics will be like that for the fer- 
romagnet with domain walls initially at every extremum 
of the potential. The decay of the magnetization for log- 
times up to Tj is given by the difference between the 
probabilities that a spin has flipped an even or an odd 
number of times. We obtain (Si(t)) ~ i{t)~ x with A = \. 
Note that this value of A saturates the lower bound of 
d/2 in contrast to the pure ID Ising case which saturates 
the upper bound of A = d M. For the symmetric RFIM 
one similarly finds that (Si(t)Si(t')) ~ (I(t') fl(t))i . 

We next study "persistence" properties. One must now 
carefully distinguish between the effective dynamics (i.e 
the walker jumping between valley bottoms ) and the 
real dynamics. The probability N(t) that a single walker 
has never crossed its starting point x(0) = xq between 
and t can be found by placing an absorbing bound- 
ary at Xo and using methods similar to the calculation of 
the endpoint magnetization in the RTFIC [p^[ . We find 
N(t) ~ l(t)- 01 at large times with X = | (c.f. X = 1 
in the pure case). A related quantity, M(t), is the frac- 
tion of starting points, Xq, for which the thermally av- 
eraged position (x(t)\x(0) = Xq) never crosses Xq up to 
time t. While in a single "run" in a given environment 
the walker typically crosses its starting point many times 
while trapped in a valley, averaging over many runs in 
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the same environment yields a (x(t)) which crosses Xq 
exactly once each time the bond on which xq lies is dec- 
imated since this causes its valley bottom to cross Xq- 
At long times, the probability M(t) thus reflects the ef- 
fective dynamics, in particular the probability that the 
bond on which xq lies has never been decimated before 
time t yielding M(t) ~ l(t)- 01 with 9 X = $=^. Indeed, 
in the biased case, the probability of no return of (x(t)} 
is like the spontaneous magnetization in the RTFIC, i.e. 
M(t) ~ \Sf for small S with = 2^1 @. 

More generally, in the effective dynamics, the prob- 
ability of exactly n crossings of the origin up to time 
t scales as ln(lni) in the unbiased case. The rescaled 
variable g = n/ln(lni) has a multifractal distribution 
Prob( 5 ) ~Z(t)-^(s) with 



Remarkably, this 8(r) is very close numerically to half 
the exact pure system result 20 \8 pure (r) — — g + 
\ (arccos( .Tr 1 ; ) ) 2 . 

To conclude, we have applied a RSRG method to ID 
random walks in the presence of static random forces and 
obtained exact results for coarsening dynamics, diffusion 
reaction models, and aging phenomena; surprisingly by 
simpler means than for the corresponding pure models. 
Extensions and details of the present results will be given 
in g7|. 

This work has been supported in part (DSF) by the 
NSF via grants DMR 9630064, DMS 9304580 and Har- 
vard University's MRSEC. 
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Since #1(5) = 0, g = | with probability 1 at large times. 
For a given walker, S(t) = j f* x(t)cLt will typically be- 
have like (x(t)). We conjecture that the probability of 
n = gln(lni) sign changes of S up to time t decays with 
the same exponent 81(g) for g < i. For larger g, the 
behavior is dominated by rare valleys with almost de- 
generate minima on opposite sides of the origin which 
yield extra sign changes in S(t). 

The persistence properties of the RFIM can similarly 
be analyzed. The probability II(t) ~ Z(i)" 6 ' that a given 
spin at x = has never flipped up to time t is equal to 
the probability that neither the nearest domain wall on 
one side, nor the nearest (opposite type) domain wall on 
the other side have crossed x — 0. Assuming the na- 
ture of the asymptotic state is as described earlier, we 
find 9 = 26>i = 1 (c.f. 9 = § in the pure case @). 
In contrast, the decay of the probability that an ini- 
tial domain survives up to time t, is S(t) ~ l(t)~^ with 
^ = 3^/5 _ 0.191 (c.f. V = 0.252 in the pure case pi 



Finally, one can study a broad class of reaction diffu- 
sion models where all particles diffuse in the same un- 
biased energy landscape and react or annihilate upon 
meeting, for example, identical particles A which react as 
A+A^A with probability 1 — r or annihilate A + A — > 
with probability r. The fraction of the valleys with no 
particle in them tends to p$, the stationary probability 
for the reaction process upon merging two valleys. Gen- 
eralizing the absorbing boundary method we obtain that, 
very generally, the probability that x — has not been 
crossed by any particle up to t decays with the exponent 
9 = 1— pq, in our example, 8(r) = 1/(1 + r). It is in- 
teresting to note that the corresponding exponent 9{r) 
associated with the thermally averaged particle trajecto- 
ries is the solution of the hypergeometric equation |i~7j| : 

9U(-r/(l + r), 28, 1) = U(-r/(l + r), 28 + 1, 1). (6) 
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